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A general research on chiral symmetry restoring phase transitions at zero temperature and finite 
chemical potentials under electrical neutrality condition has been conducted in a Nambu-Jona- 
Lasinio model to describe two-flavor normal quark matter. Depending on that mo/A, the ratio 
of dynamical quark mass in vacuum and the 3D momentum cutoff in the loop integrals, is less or 
greater than 0.413, the phase transition will be second or first order. A complete phase diagram of 
u quark chemical potential versus mo is given. With the electrical neutrality constraint, the region 
where second order phase transition happens will be wider than the one without electrical neutrality 
limitation. The results also show that, for the value of mo/A from QCD phenomenology, the phase 
transition must be first order. 
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I. INTRODUCTION 



Nambu-Jona-Lasinio (NJL) model [l| is a good low 
energy phenomenological model which can be used to 
simulate Quantum Chromodynamics (QCD) and is suit- 
able for researching spontaneous breaking of symmetry 
and its restoring at finite temperature and finite chemi- 
cal potential in quark matter [3, I n a realistic research 
of the normal quark matter without color superconduc- 
tivity, owing to that the quarks carry electrical charges, 
one must impose electrical neutrality condition on the 
quark matter. In this paper, we will generally examine 
symmetry restoring phase transitions at zero tempera- 
ture and high density in a NJL model with two-flavor 
normal quark matter under electrical neutrality condi- 
tion, however, the confinement problem which is beyond 
the power of a NJL model will not be touched on. This 
examination will make us get more deep-going under- 
standing of the affect of electrical neutrality on phase 
transitions and certainly has important theoretical sig- 
nificance. It is noted that such general research has not 
yet appeared in present literature. The preceding rel- 
evant researches either had completely not touched on 
electrical neutrality requirement [i[ or only limited rel- 
evant parameters e.g. the momentum cutoff A and the 
dynamical quark mass mo in vacuum to some given val- 
ues by phenomenology. § . 

It is a well-known fact that in the vacuum of QCDjChi- 
ral symmetry is inevitably spontaneously broken [l|, 0, 0, 
IU, so in our model, we will first assume that the quarks 
in vacuum have acquired a common non-zero dynami- 
cal mass from the quark-antiquark condensates, then re- 
search that, as increase of the quark chemical potential, 



the dynamical quark mass could change to zero and the 
chiral symmetry will finally be restored . In the process, 
the electrical neutrality condition will always be main- 
tained and the order of the symmetry restoring phase 
transitions will be the key point of our concern. 

For keeping electrical neutrality of the quark mat- 
ter, we will add the contribution of free electron gas to 
the effective potential. The discussions will be made in 
the mean-field approximation. Throughout the paper, a 
three-dimension momentum cutoff will be used. 

The paper is arranged as follows. In Sect[TT]we will give 
the effective potential of a NJL model describing two- 
flavor normal quark matter, its extreme value equation 
and the electrical neutrality condition. In Sect lIII| we 
will discuss the second order phase transition in small dy- 
namical quark mass in vacuum, or equivalently, in large 
momentum cutoff, and in Sect HVI expound the second 
and the first order phase transition in larger ratio of the 
above two parameters. Finally, in Sect|V]we come to our 
conclusions. 



II. EFFECTIVE POTENTIAL, GAP EQUATION 
AND ELECTRICAL NEUTRALITY CONDITION 

The Lagrangian of the NJL model describing the two- 
flavor normal quark matter can be expressed by 



C = qi^d^q + G s [{qq) 2 + {qi^rq) 2 



(1) 



with the quark Dirac fields q in the ST// (2) doublet and 
the SU C {2>) triplets, i.e. 



i = r,g,b, 
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where the subscripts i — r,g,b denote the three col- 
ors (red, green and blue) of quarks, r = (ti,T2,Ts) 
are the Pauli matrices, Gs is the four-fermion coupling 
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constants and we have omitted the bare mass of the 
quarks. Assume that the four-fermion interactions can 
lead to the scalar quark-antiquark condensates (qq) = <f), 
then the chiral iSC//x(2) <g> SUfn(2) flavor symmetry of 
the Lagrangian (1) will be spontaneously broken down 
to SUfv{2) and the quarks will get dynamical mass 
to = —2Gs(qq) ■ In the mean field approximation 
we can write the effective potential of the model in the 
T — > limit by 

V(m,/i,Me) = 

+ - £ P )(M« - E p ) + (/!„ -> pa)]} - 77TI (2) 



where E p = yjp 2 + to 2 , fi = —dV/dn is the quark chem- 
ical potential corresponding to the total quark number 
density n, [i e is the chemical potential of electron and 

2 1 

/in = /i - -fie, fJ-d = V + gA*e = to + to (3) 

are respectively the chemical potentials of the u and d 
quarks. The second equality in Eq.(3) is usually refers 



as beta equilibrium [9]. In vacuum we have fx = fx e = 0, 
thus the effective potential (2) is reduced to 

™-iB;- a J 

With a 3D momentum cutoff A, we may find out that the 
extreme value points of Vo are 1) m = and 2) to = toq, 
where Too obeys the gap equation 

(4) 

ft is easy to verify that if 

1/2G S < 3A 2 /tt 2 , (5) 

then m — will be a maximum point and simultaneously 
Eq.(4) will have non-zero solution too which is a mini- 
mum point. This means spontaneous breaking of chiral 
symmetry in vacuum. This will assumedly be our pre- 
supposition of discussions in this paper. Hence we may 
replace l/2Gg in the effective potential V(m, fi, to) by 
using Eq.(4). As a result, we obtain from Eq.(2) 



V(m,fj,,fi. 

3 
4^2 



to 2 I 2A A /A 2 + to 2 , - 2TOgln 



A+ ^A 2 +ml 



to 



+ 2A 4 - A \/ A 2 +to 2 (2A 2 + to 2 ) + to 4 In 



4, A + VA 1 ^ 



-4^2 i - w ) 



77? 



f27T 2 



(6) 



For deriving electrical neutrality condition, it is noted 
that the electrical charge density in the two-flavor quark 
matter with electrons is 



f 



from which we may obtain 



-n d 



to 



dV 

dn e 



dV dn Q 
driQ dn e 



Hence the electrical neutrality condition will become 
tiq = —dV/dfiQ = dV/dfie = and has the following 
explicit expression 



dV 
dfj, e 



1 

3^ 



29(ij u ~ m)(fi 2 u -to 2 ) 3 / 2 



-6(to~m)(p 2 d -rn 2 ) 3/2 -d 



0, 



to = Vu + to- (J) 

Eq.(7) is a restraint condition about to, fi u and /i e . In 
fact, in the effective potential (6), instead of \i and /u e , 
we can first consider fi u and fj, e as two chemical potential 
variables, then by Eq.(7), only one of them, the u-quark's 
chemical potential is left as a single independent one, 
since \i e may be viewed as a function of ji u and to by 
Eq.(7). This treatment will bring about great conve- 
nience for the discussions of phase transitions. 

For research in ground state of the system, we must 
consider the extreme value points of V(m, fx U i Me) deter- 
mined by the equation 
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dV i_ 

dm 7T 2 



Ai/A 2 + ml - ml In A + v ' ^ + ^ - (to,, - - m ) 



m 



/)/ N / /~~2 2 2, Mu + VMS - ™ 2 

0(/x„ - m) I //„ V M« -m 2 -m In ^— 



+ (m« -» Md) 



and second derivation of V(m, /U e ) over m under the constraint given by Eq.(7) 



d 2 V d 2 V ( d 2 V 



dm 2 dm 2 \dmd/j, e ) dfi 2 



,d 2 V 



with 



8 2 V 
dm 2 



= \ { A JXM^ - ml In A + ^ +m K 3m 2 ( In 

7T 2 V mo ^ 



A + VA 2 + to 2 A 



+ - 



0(/z u - m) I fi u fi 2 u -m 2 + 3m 2 In 



+ (m« -» Md) 



(8) 



(9) 



A 3 



VA 2 + m 2 I VA 2 + m 2 



(10) 



and d 2 V / dmdfi e and d 2 V/dfi 2 can be obtained from 
dV/dfi e in Eq.(7) when Eq.(3) is taken into account. 

Extremal feature of the point m = is quite important 
for determination of behavior of the effective potential V. 
Substituting m = 0, which is obviously a solution of the 
extreme value equation (8), into Eq.(7), we will obtain 
the electrical neutrality condition at m = 



0. 



(11) 



Let ?7 = /i e //i M , then from Eq.(ll) we may obtain a real 
number solution rj = 0.256. Furthermore, from Eqs.(9) 
and (10), second derivative of V over m at m — be- 
comes 



d 2 V 



dm 2 



where 



d 2 v 



m=0 



dm 2 



m=0 



3 
2^ 2 



£-[1 + (1 + ^)2] (^.^ 



(12) 



2 

M«c — 



2 / 2 A+y/F + 7% 

i+(i +?7 ) 2 r° ln - ^ 



mo 



-A^/a 2 + to 2 + A 2 ^) . 



(13) 



Hence, m = will be a maximum (minimum) point of 
V if /i 2 < /i 2 c (/x 2 > /i 2 c ), and the extremal feature of 
m = will be determined by derivatives of higher order 
of V over m if fi 2 u = /x 2 c . 



III. SECOND ORDER PHASE TRANSITION IN 
SMALL mo /A 

It is easy to verify that when fi u — \i e — (or equiv- 
alcntly, fi = fi e = 0), the effective potential given by 



Eq.(6) will reproduce the spontaneous chiral symmetry 
breaking in vacuum and the quarks get the dynamical 
mass to = to - Starting from this, we will first analysis 
possible chiral symmetry restoring by second order phase 
transition as the chemical potential fi u increases and de- 
rive the fi u — too critical curve of the phase transition. 
The variation of V(m, fi Ul fi e ) as increase of fi u will be 
discussed successively. 

1) < fi u < too < fi U c- Here we impose the limitation 
too < fi uc which, by Eq.(13), implies that A/too > 2.865, 
i.e. we are confined to the region too/A < 0.349. In 
this case, from Eqs.(7)-(12), the effective potential V 
will have the minimum point to = too with fi e = 
which comes from the electrical neutrality condition for 
fi u < to, and the maximum point to = with the 
electrical neutrality constraint fi e = rjfi u for fi u > to. 
A following question is that, for the case of fi u > to, 
whether there is any other non-zero solution satisfying 
the electrical neutrality condition (7) for the gap equa- 
tion {dV/dm)/m = coming from Eq.(8) ? To answer 
this question, we note that in this case the above two 
equations may be reduced to 



2(/4 - to ) 



2\3/2 



(Ma 



TO 2 ) 3 / 2 = £ 



(14) 



and 



TOq In 



A + v/A 2 + ml 



toq 



1 



to In • 



(A + VA 2 + to 2 ) 2 



2 {fiu + V 74 - m 2 )(fi d + ^Jfi 2 d - to 2 ) 

+ ^(M«V / M« - TO 2 + fj, d yl 1 fl 2 - TO 2 ) - A^A 2 +TO 2 . 

(15) 



4 



Denote that 

a = A/mo, x = fi u /mo, /? = m//i u , a = // e /A*« 
then Eqs.(14) and (15) can be changed into 

2(1 _ ^2)3/2 _ [(1+a) 2_ /3 2 ] 3/2 =a 3 

and 

x 2 /3 2 



(16) 



ln(a + + 1) - avV + 1 



In- 



(a + v/a 2 + x 2 f3 2 f 



2 a; 2 (l + v^)[l + a + v/(l + a) 2 - /3 2 ] 
+ y [v 7 !^ 2 + (1 + a)v/(l + a) 2 -/3 2 ] - a^a 2 + x 2 (3 2 



(17) 



r 



It is easy to check that when a > 2.865, Eqs.(16) and 
(17) have no solution with x < 1 and (3 < 1. In other 
words, when a -1 = too/A < 0.349, the gap equation 
(dV/dm)/m = and the electrical neutrality equation 
(7) have no solution with fi u < to and m < fi u indeed. 

To sum up, when fi u < to < fi UCl the effective po- 
tential will have only a maximum point to — and a 
minimum point m — ttiq accompanied with fi e = 0. The 
latter corresponds to ground state of the system which is 
similar to the case of vacuum. This shows that sponta- 
neous chiral symmetry breaking in vacuum will be main- 
tained in the case of fi u < too and \i e = 0. The dynamical 
quark mass m could be changed only if fi u > m . 
2) mo < fi u < fiuc- In this case, m = is still a maxi- 
mum point of V by Eq.(12). On the other hand, it is not 
difficult to see that now Eqs.(7) and (8) have no solution 



with fi u < to. Thus we are left only the case of fi u > to 
and the gap equation (dV/dm) /to — and the electrical 
neutrality condition will take the forms of Eqs.(14) and 
(15). In view of the definition of fi uc given by Eq.(13) and 
the constraint /i 2 < fi uc , it may be deduced that Eqs.(14) 
and (15) have non-zero solution (mi,/j e i). From Eq.(9) 
we obtain that 



d 2 V 



dm 2 



d 2 v 



(roi,|Uel) 



dm 2 



( d 2 V 
\dmdfi e 



,d 2 V 



dp 2 



(™l^el) V ~r~e (mii(Ue 

(18) 

The second term in the right-handed side of Eq.(18) is 
positive since d 2 V/dfi 2 derived by Eq.(7) is always neg- 
ative. The first term may be calculated with the result 



d 2 v 

dm 2 



(A + v/A 2 + m 2 ) 2 



(mi.Mei) 



— I —L In 

^ 1 2 11 (Pu + y/nl - m 2 )(fi d + ^fi 2 - m{) A 2 + m{ 



to 2 . In 



to 2 A 



A + y/A 2 + to 2 



too 



A^ A 2 + to 2 + A^ A 2 + m\ - hjj^yjul - m\+ fi ds ] fi\- m\) 



sjh 2 + mj 

> ^2 | + A) - \(Vu\jfi 2 a - m\ + fi d ^fi 2 - to 2 ) - A 2 ^1 



0V 2 + to 2 



>0, 



Mel 



where we have used Eq.(15) with to replaced by m\ and In this way, it is proven that 
the condition /z 2 < fi uc , or equivalently, 



i[l + (l + ?7 ) 2 ] M 2 c > l -[l + {l + rff]nl > + 



Mel 



d 2 V/dm 2 \ (mi ^ el) > 0. 

Since for a fixed fi u , {m\,fi e i) is now the only minimum 
point of V, so it will correspond to the ground state of 
the system satisfying electrical neutrality. 

It may be found by examining Eq.(15) that when fi u — 



■5 



mo and /x e = 0, we have m — jtiq. As fi u increases from 
Too and when Too < \i u < \x uc and fi e > 0, the left-handed 
side of Eq.(15) keeps unchanged, and in the right-handed 
side, the second and the third term will change from zero 
to positive, then the first term has to decrease so that 
the original Too will change to to < too, and going up 
of the fourth term — A\/A 2 + m 2 is consistent with the 
reduction of to in the first term. This means that to will 



decrease from Too and finally it may continuously reduce 
to zero, thus we come to a critical point of second order 
phase transition at which the broken chiral symmetry will 
be restored. The second order \i u — m critical curve will 
be denoted by C2(to ) whose equation can be obtained 
by setting to = in Eqs.(14) and (15) and has the explicit 
expression 



l + (l + r/) 2 



In 



too 



+ TO, 



o'+A 2 



1/2 



C 2 (mo). 



(19) 



3) mo < IJ, U = fi U c- In this case, the extreme value equa- 
tion (8) has the only solution to = and the electrical 
neutrality condition (7) gives \i e — rjfi u with r\ — 0.256. 
At the only extreme point to = of V, the n-th deriva- 
tive of V over to under the electrical neutrality condition 
expressed by Eq.(14) may be found out to be 



0, 

9 In A 2 

M„ C G(?7) 



when n = 2, 3, 5 
when n = 4 



m=0 



15 / 9 



^•ffi^ ) , when n = 6 



where we have used the denotation 

G(r ? ) = (l+,)exp ^ + 2 ^ + v 

Eq.(20) implies that when 

ML< A 2 G- 1 ( ?? ), 



(20) 



(21) 



to = will be the only minimum point of V and the 
broken chiral symmetry will be restored in a second order 



phase transition. In view of Eq.(19), the condition (21) 
also means the constraint 



toq 
A 



< 0.342, 



(22) 



i.e. in the region < too/A < 0.342, /i u — /i uc is a 
curve of second order phase transition. It is noted that 
the constraint (22) is consistent with the presupposition 
mo < Hue, i-e. too/A < 0.349. 

4) fi u > ii uc . We still confine ourselves to the case of 
a4c < A 2 G" 1 (r7). Obviously, by Eq.(12) and fi u > fj, uc , 
to = is now a minimum point of V. But one can 
raise such a question that in the above condition, whether 
dVjd\x e = and (dV/dm)/m = also have some non- 
zero to solutions with m < /j, u ? To answer this query, we 
can rewrite the above equations, i.e. Eqs.(14) and (15) 
by 



2(1 _ 02)3/2 _ [(1 + a) 2 _ 02J3/2 



(23) 



and 



I[l + (l + ?7) 2 ] = ^!ln - (^V^V) 2 

2 L 2 7 2 (1 + y/T=p){\ + a + + a)* - /3 2 ) 



+ 



7 



VT? + (1 + a)y/(l + a) 2 - 0* +6 2 -6 v /& 2 ~fW, 



where we have used the denotations 



a = fj, e /fi u , = m/fi, u , 7 = Mu/Muc 

b 2 = A 2 //i 2 c = [1 + (1 + V f] a 2 /2 \ln(a + Va 2 + 1) - a^a 2 + 1 + a 2 



A/m . 



(24) 



(25) 



It turns out by numerical calculation that Eqs.(23) and (24) may have solutions with [3 < 1 (to < n u ) and 
7 > 1 {ji u > ji c ) except that a — A/to is quite large (e.g. a > 5). In other words, when /i > /j, uc , the effective 
potential V could have extreme value point with m ^ 0. However, it can be proven that to = is always the least 
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FIG. 1: The fj,„ — mo phase diagram of the model. Here [i u and mo are both scaled by the momentum cutoff A. The straight 
line through the origin and the point C is ^ u = mo. The critical curve ^ u = C*2(mo) of second order phase transitions starts 
from the origin, through the points A and C, ends at the point B. The critical curve fi u = Ci(mo) of first order phase transition 
begins from the point B then extends to the right in the region above jj, = flua- The B is a tricritical point. 



minimum point. In fact, by a direct calculation, we may obtain the difference between the values of V at the extreme 
value points m/0 and m = expressed by 



V\ m ^ - V\ 



rn—0 



3 

4^ 



12tt 



> g§<3 



|2A 4 - 2AVA 2 +rn 2 + A 2 m 2 - |[1 + (1 + r,) 2 ] 
y {3/x„(/4 - m 2 f' 2 + 3( Mu + /x e )[(M« + Me) 2 - m 2 ] 3 / 2 + M 4 } 

[3/4 + 3(M« + Me) 4 + Me] | Me=wt „ 

2b 4 - 2b^b 2 +f3 2 1 2 + 6 2 /3 2 7 2 - 1 + (1 n +??)2 /3 : 



12tt 2 
1 



-OtfJbu 



(1_ ^2)3/2 

(1 + a)[(l + a) 2 - /3 2 ] 3/2 - a 4 /3 + 1 + (1 + f]) 4 + ij 4 /^ , when fj, u > fi u 



r 



(26) 



In the condition ^ 2 C < A 2 G _1 or b 2 > G(r]), substituting 
all the possible solutions of a,/3 and 7 obtained from 
Eqs.(23) and (24) into Eq.(26), we will always obtain 
that T^|m^o — ^|m=o > 0. This indicates that m = 
is indeed the least minimal value point of V and when 
Mil > Muc) the chiral symmetry has been restored 
through a second order phase transition in the case of 
/i 2 c < A 2 G~ 1 . FIG. Q]is the complete fj, u — mo phase 



diagram of the model. In this diagram, the discussed 
second order phase transition above will correspond to 
the segment of the curve fi u = C^mo) from the origin to 
the point A whose location is determined by the equality 
C 2 (m ) = KG- 1 ' 2 ^). 
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IV. SECOND AND FIRST ORDER PHASE 
TRANSITIONS IN LARGER mo/A 



termine the critical curve of a first order phase transition. 
In electrical neutrality condition, these equations read 



It may be seen from Eq.(20) that, when going along the 
curve fi u = [i uc toward the region with [i} uc > A 2 G~ 1 (rj), 
or equivalently, in view of Eq.(22), mo/A > 0.342, one 
will get m = becoming a maximum point of V, how- 
ever when fi u > fj, UC} to = is again a minimum point. 
Such change of minimax property of to = could lead to 
two possibilities: either a second order phase transition 
will continue or a first order phase transition will hap- 
pen. For examining a concrete realization of the above 
two possibilities, we will start from the equations to de- 



dV 

V(m — 0) = V(m = mi), — — 

o\x e 



dV_ 

dm 



jm 



0, 



d 2 V 



dm 2 



m=0 

> 0. 



dV 

Ofl e 



0, 



(27) 



Since dV/dfx e \ m —o — determines only the ratio 77 = 
Me/A*«c at m = 0, Eqs.(27) will have the following explicit 
expressions: 



J 



/4 



[I + (1 + r?) 4 + 77 4 /3] = m\ 



-2A 4 + 2A d i/A 



+ [0(jiu ~ mi)fi u {nl - m?) 3 / 2 /3 + {fi u -> fi d )] + (i 4 e /9, 
29(n u - mM - m\fl 2 - 6([i d - mi)(^ - m\fl 2 -fi=0, 



(28) 
(29) 



l + (l + T)f 



2 



A^A 2 + ml + ml 



In 



A + y/N 2 + ml 




mi) 



mi 

[l\ — ml — ml In 
I 



m\ 



+ (fX u -> /l d ) 



(30) 



and 



2 > 2 



(31) 



We will discuss respectively the two cases of [i uc > 
toq and fi uc < mo which could appear when ^t„ c > 

1) fi u > fi uc > too. Since jj, uc > too can be satisfied 
only if toq/A < 0.349 and in view of Eqs.(21) and (22), 



the limitations fi uc > to and fi 2 c > A 2 G~ 1 (rj) will corre- 
spond to the region 0.342 < too /A < 0.349 in the fi u — too 
plane, i.e. the region between the points A and C in 
FIG.[T] In the present case of y, u > Too, the gap equation 
(30) has no solution for too < /i u < mi, so we need to 
consider only the solution for fi u > mi, thus all the 6- 
functions in Eqs.(28)-(30) may be removed. As a result, 
Eq.(28) will be changed into 



-[l + (I + 7 ? ) 4 + r, 4 /3] = /3V 



1 + (1 + r/) 2 



-b 2 



2& 4 + 2&V& 2 + /3 2 7 2 



+y {(1 - /3 2 ) 3/2 + (1 + «)[(! + a) 2 - /3 2 ] 3/2 + a 4 /3} 



(32) 



r 



Eqs.(29) and (30) will be separately identical to Eqs.(23) 
and (24), and Eq.(31) will simply become 7 > 1, where 
we have again used the denotations given by Eq.(25). By 
numerical solution of Eqs. (23), (24) and (32), we have 
found that in the whole region 0.342 < toq/A < 0.349, 



there is no solution with toi 7^ 0. This means that in this 
region a first order phase transition could not happen. 
In fact, the obtained solution is toi = 0, which makes 
Eq.(32) becomes a trivial identity and Eqs. (23) and (24) 
are reduced to fJ, e /fJ>u = '7 an d A* = AW- The last two 
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equalities are precisely the equation of the curve C^mo). 
Hence we can conclude that in this region, the phase 
transition remains to be second order. 
2) Hue < M« < m o- In FIG. [1] this corresponds to the 
right-handed side of the point C. For m < fi u , it is 
easy to verify that the electrical neutrality equation (23) 
and the gap equation (24) have non-zero solution mi for 



Mit > Muc and n uc > A 2 G~ 1 (rj) and m = is a minimum 
point. This leads to the possibility to generate a first 
order phase transition. For m > fi u , by solving Eqs.(29) 
and (30), it is obtained that m = mo with fi e = is 
a minimum point. Hence, with mi — mo and fi e = 
being taken, Eq.(28) which determines first order phase 
transition point may be changed into 



Mu = 



l + (l+77) 4 +7? 4 /3 



mi In 



A+ y/A^T 

m 



mnA 



A 2 



2A 4 + 2A 3 A /A 2 



1/4 



Ci(m ). (33) 



Eq.(33) expresses the equation of the curve Ci(mo). 
Ci(rafl) may becomes a first order phase transition curve 
only if /i u = Ci(mo) > [i U c- From this constraint we 
obtain 

m /A > 0.413 (34) 

which is obviously in the right-handed side of the point 
C. At mo/A = 0.413, the curve /i u = Ci(mo) and the 
curve fi u = C*2(mo) intersects. The intersection point 
denoted by B becomes the starting point of the curve 
M« = Ci(m ) above fi > [i uc (m ). 

In the region 0.349 < mo/A < 0.413, i.e. in the seg- 
ment between the points C and B, similar to the case of 
fi uc > m , it can be proven that Eqs.(23), (24) and (32) 
have only the solution mi = 0, Me/Mu = V an d M = M«c, 
i.e. the solutions are actually reduced to the second order 
phase transition curve fi u — C2(mo). Therefore, in the 
C — B segment, we still have second order phase transi- 
tion represented by the critical curve fi u = C2(mo). 

In summary, in the region with /i 2 c > A 2 G~ 1 (ry) or 
mo/A > 0.342, the critical curve /i u — /i uc = C2(mo) 
of second order phase transition may be extended to 
the point B where mo/A = 0.413, then a critical curve 
[i u = Ci(mo) of first order phase transition will start 
from the point B in the region of fi u > fj, uc . So the point 
B is a tricritical point. 

It should be indicated that in the case without and 
with electrical neutrality constraint, the feature of phase 
transition of the NJL model is different. Without electri- 
cal neutrality condition, as was discussed in Ref.[4], the 
second order phase transition curve will end at a similar 
point A, then from A through the point C straight to the 
right-handed side of the point B, one will always have a 
first order phase transition curve. However, in present 
case with electrical neutrality requirement, from A to B, 
one continue to get a second order phase transition, in- 
stead of a first order one. 

The total conclusions of this paper come from a gen- 
eral research of the used NJL model, i.e. the parameters 



of the model, the 3D momentum cutoff A and the dy- 
namical quark mass mo in vacuum (correspondingly, the 
four-fermion coupling constant Gs) have been considered 
as arbitrary ones. If the discussed NJL model is used to 
simulate QCD for normal quark matter and the conven- 
tional phenomenological values of the parameters Q 

G s = 5.0l63GeV~ 2 ,A = 0.6533G e y 

are taken, then by Eq.(4) we will obtain mo/A = 0.48. 
Based on the results of present paper, chiral symmetry 
restoring at high density in this model must be first order 
phase transition. 



V. CONCLUSIONS 

In this paper, we have generally analyzed chiral sym- 
metry restoring phase transitions at zero temperature 
and high density in a NJL model to describe two-flavor 
normal quark matter under electrical neutrality con- 
dition. It has been found that the feature of phase 
transitions is decided by the ratio mo/A, where mo is 
the dynamical quark mass in vacuum and A is the 3D 
momentum cutoff of the loop integrals. Depending on 
mo/A is less or greater than 0.413, the phase transition 
will be second or first order. As a comparison, the 
resulting region where a second order phase transition 
happens is wider than the one in the case without 
electrical neutrality constraint. For the value of mo /A 
based on QCD phenomenology, the phase transition 
must be first order. 

The present discussions of normal quark matter based 
on a NJL model can be generalized to the case of color 
superconducting quark matter where one must also con- 
sider diquark condensates, besides the quark- antiquark 
condensates, and at the same time, also impose color 
neutrality condition, besides electrical neutrality one. 
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